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Table 1. Basic parameters of NGC 4753. 



Position: a(2000) = I2 h 52 m 23.23 s ; (5(2000) = -01 12 00.50 
Distance = 8.7 Mpc. 
Velocity = 1225 km s~ x 

Magnitude : U T = H-26 ± 0.10; B T = 10.85 ± 0.10; V(Johnson) = 10.30. 
FIR flux densities : /(12/xm) = 340 ± 42 mJy; /(25/ito) = 310 ± 71 mJy; 

/(60/im) = 2640 ± 60 mJy; /(100/wn) = 8010 ± 176 mJy. 
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Table 2. Transformation Coefficients. 



Quantity 


a 


/3 


V-Vq 


0.399 ± 0.004 


0.021 ±0.006 


B-V 


-0.473 ±0.015 


1.416 ±0.018 


V -R 


-0.041 ±0.008 


1.001 ±0.017 


B - R 


-0.482 ± 0.008 


1.248 ± 0.006 


Table 3. 


BCES and OLS linear rej 


jression coefficients. 


fit 


slope 


intercept 


BCES(A B |A y ) 


2.61 ±0.476 


-0.292 ±0.102 


OhS{A B \A v ) 


1.20 ±0.046 


0.009 ±0.011 


BCES(A y |A s ) 


0.747 ±0.113 


0.106 ±0.0255 


OLS(Ay|A B ) 


1.323 ±0.054 


-0.017 ± 0.013 


BCES bisector 


1.32 ±0.0881 


-0.0168 ±0.0197 


OLS bisector 


1.26 ±0.048 


-0.003 ±0.012 


VCVS{A R \A V ) 


2.10 ±0.493 


-0.315 ±0.110 


BCES(Av|A fl ) 


0.188 ±0.131 


0.112 ±0.0294 


BCES bisector 


0.77 ±0.0811 


-0.0182 ±0.0184 


BCES(A B |A fl ) 


4.53 ±2.08 


-0.404 ±0.314 


BCES(Ar|Ab) 


0.852 ±0.146 


0.0221 ±0.0254 


BCES bisector 


1.66 ±0.171 


0.0221 ±0.0254 
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Table 4. Parameters used in the estimation of dust mass. 



Parameter Value References 



AsUicate lO" 2511 cm 2 - 5 /H 1 

AgrapMte IQ- 2516 Cm 2 5 / H 1 

n H 20 cm- 3 2 

Pd,graphite 3.26g CTU^ 3 1 

Pd.silicate 3.3g cm^ 3 1 



References for Table 4. 
1. Draine & Lee (1984) 2. Hoyle & Wickramasinghe (1991) 



Table 5. Masses of different components of ISM in 
NGC 4753. 



Component Mass (Mq) 

a DuSt {M d . op Ucal) 1.5 x 10 5 

°Dust (M d . /JM s) 3.46 x 10 5 

fe HI < 5.28 x 10 8 

a Molecular gas 0.2 - 2.4 x 10 8 

c Ha gas 0.67 x 10 5 



a derived value 

b based on 21 cm line emission observed by Knapp et al. 1985 
c taken from Singh et al. 1995 
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Dust Properties of NGC 4753 
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ABSTRACT 

We report BVR surface photometry of a lenticular galaxy, NGC 4753 with promi- 
nent dust lanes. We have used the multicolor broadband photometry to study dust- 
extinction as a function of wavelength and derived the extinction curve. We find the 
extinction curve of NGC 4753 to be similar to the Galactic extinction curve in the vis- 
ible region which implies that the sizes of dust grains responsible for optical extinction 
are similar to those in our Galaxy. We derive dust mass from optical extinction as well 
as from the far infrared fluxes observed with J/M^Q. The ratio of the two dust masses, 
iH d ' IRAS , is 2.28 for NGC 4753, which is significantly lower than the value of 8.4 ± 1.3 

^"d, optical 

found previously for a large sample of elliptical galaxies. The total mass of the observed 
dust within NGC 4753 is about a factor of 10 higher than the mass of dust expected 
from loss of mass from red giant stars and destruction by sputtering and grain-grain 
collisions in low velocity shocks, and sputtering in supernova driven blast waves. We 
find evidence for the coexistence of dust and Ha emitting gas within NGC 4753. The 
current star formation rate of NGC 4753, averaged over past 2 x 10 6 yr, is estimated 
to be less than O^IMqi/t' 1 . A substantial amount of dust within NGC 4753 exists in 
the form of cirrus. 



Subject headings: galaxies: elliptical and lenticular, cD — galaxies: ISM — galaxies: 
photometry — infrared: galaxies — (ISM:) dust, extinction 
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1. Introduction 

It has now been well established that the 
presence of dust in early-type galaxies is the 
rule rather than an exception. Cool interstel- 
lar dust has been observed by IRAS (Jura 1986) 
and by means of optical extinction (Hawarden 
et al. 1981; Sadler & Gerhard 1985; Ebneter et 
al. 1988; Goudfrooij et al. 1994b). Using IRAS 
data, Knapp et al. (1989) have shown that about 
45% of the ellipticals and 68% of the SO galaxies 
have been detected at 60 /im and 100 //m bands 
while Goudfrooij et al. (1995) have reported a de- 
tection rate of 61% for their sample of Shapley- 
Ames elliptical galaxies. The Hubble Space Tele- 
scope (HST) survey of ellipticals in the Virgo 
cluster has revealed the presence of dust in the 
nuclei of almost every galaxy (Jaffe et al. 1994). 

Apart from dust, early-type galaxies also con- 
tain substantial amounts of gas at various tem- 
peratures. Cold gas (T < 100 K) has been ob- 
served via its HI 21 cm line emission (Knapp 
et al. 1985) and CO line emission (Wiklind k, 
Rydbeck 1986; Phillips et.al 1986); warm gas 
(T rs 10 4 K) via optical emission lines (e. g. 
[Oil], [OIII]) (Caldwell 1984; Phillips et al. 1986) 
and hot gas via their X-ray emission (Forman 
et al. 1985; Canizares et al. 1987; Fabbiano et 
al. 1992). The interstellar medium (ISM) of 
early-type galaxies, thus contains many differ- 
ent phases from cold to very hot. Typically 
they contain an amount of cool dust of mass 
~ 10 4 - 1O 6 M (Jura et al. 1987; Goudfrooij et 
al. 1994c), hot gas of mass ~ 10 9 - 1O 1O M (For- 
man & Jones 1985; Fabian et al. 1987), and small 
amounts of warm gas of mass ~ 10 3 — 10 5 Mq 
(Phillips et al. 1986; Caldwell 1984). 

The study of properties of different phases of 
ISM in early-type galaxies provides important 
clues to the origin, nature and fate of the in- 
terstellar matter. In particular, dusty early-type 
galaxies provide suitable environment to study 
the nature of extra-galactic dust grains, either 
via their far infrared (FIR) emission or optical 



extinction. The data can be used to study the 
physical mechanism operating on them which can 
determine their size distribution, temperature, 
abundances, formation, destruction etc. The 
physical properties of the dust grains are func- 
tions of time and can be used as an indicator for 
the time elapsed since the dust was last substan- 
tially replenished (see Goudfrooij et al. 1994c). 

The emission line regions in X-ray bright early- 
type galaxies are expected to be dust-free in view 
of very short life time (~ 10 6 — 10 7 yr) of dust 
grains. However, Goudfrooij et al. (1994c) found 
that the emission line regions in these galaxies 
are essentially always associated with substan- 
tial dust-extinction. This dilemma can be re- 
solved in the evaporative flow scenario (de Jong 
et al. 1990) in which the observed dust immersed 
in hot, high pressure plasma may be replenished 
by the evaporation of cool clouds. The presence 
of cool clouds is assumed to be due to the cap- 
ture of a companion dwarf galaxy, rich in dust 
and gas. A recent study of nuclear dust in 64 
elliptical galaxies imaged with HST has shown 
that the dust and gas are generally dynamically 
decoupled (van Dokkum & Franx 1995). More- 
over, the kinematics of stars and gas in ellipticals 
with large dust lanes are found to be decoupled 
(e.g. Bertola et al. 1988). These conclusions gen- 
erally indicate an external origin of dust in these 
galaxies. Further, the dynamical state of dust 
and gas can be used to probe the intrinsic shape 
of the underlying galaxy. The dust and gas, set- 
tled in the galaxy potential, are allowed in sta- 
ble closed orbits which indicates a plane in the 
galaxy (Gunn 1979; Merrit & de Zeeuw 1983; 
Habe & Ikewchi 1985, 1988). 

In this paper, we present a detailed study of 
dust properties of an early-type galaxy NGC 4753 
with very prominent dust lanes. The galaxy 
NGC 4753 is located at a distance of 8.7 Mpc 
as determined from a distance modulus of 29.7 
by Buta et al. (1985) using published estimates 
and the light curve of supernova SN 1983g ob- 
served within the galaxy. NGC 4753 has been 
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classified ELS Si peculiar SO type galaxy in the 
Hubble atlas (Sandage 1961) because its under- 
lying luminosity distribution resembles an SO 
type galaxy. However, the third reference cat- 
alogue (RC3) (de Vaucouleurs et al. 1991) clas- 
sifies it as an irregular type, 10. The complex 
dust lanes in NGC 4753 pass through its center. 
Steiman-Cameron et al. (1992) have shown using 
R band photometry that the dust lanes lie in a 
disk which is strongly twisted by differential pre- 
cession. We have carried out BVR surface pho- 
tometry of NGC 4753 and derived its color and 
extinction maps. Based on these, we estimate 
the dust mass from optical extinction. We have 
also derived the dust mass based on FIR emis- 
sion. The observed optical (RC3) and infrared 
(Knapp et al. 1989) properties of NGC 4753 are 
summarized in Table 1. The paper is organized 
as follows. In the next section, we provide details 
of our observation. In §3 we present our method 
of analysis and the results obtained. In §4 we 
discuss our results followed by conclusions in §5. 

2. Observations 

NGC 4753 was observed with the Vainu Bappu 
Telescope (VBT) on the night of March 13, 
1996. The observations were carried with a liq- 
uid Nitrogen-cooled TEK 1024 x 1024 CCD chip 
placed at the prime focus of the 2.3 m reflec- 
tor. The nominal value of f-ratio is f/3.25. The 
pixel size of 24 /im square per pixel of the CCD 
chip gives a scale of 0.6" per pixel and a total 
field of 10.24' x 10.24'. The observations were 
carried out under photometric conditions. The 
seeing (FWHM) was in the range 2.1" — 2.6". 
The exposure times were chosen not to satu- 
rate the CCD pixels because of the presence of 
bright stars in the neighbourhood of the galaxy 
NGC 4753. In order to achieve a good signal-to- 
noise ratio, two images of NGC 4753 were taken 
in each of the broad band filters B,V and R with 
exposure times 300, 180 and 60s respectively. To 
correct for the bias level, bias images were taken 
just before and after the galaxy observation apart 



from many bias frames taken during the observ- 
ing night. Several flat field images in each filter 
were taken by exposing to the twilight and dawn 
sky in order to correct for the non-uniformity of 
response of pixels of the CCD. The standard star 
field in the "dipper asterism" region of the open 
cluster M 67 was observed for photometric cali- 
bration. Both NGC 4753 and M 67 were observed 
close to the zenith. 

3. Analysis and Results 

We have used IRAF0 software package for the 
basic reduction and analysis of CCD images. Ex- 
amination of the bias frames showed that the 
mean bias level did not vary significantly over 
the observation run. An average bias frame, 
constructed from bias frames taken very close 
to the galaxy observation, was subtracted from 
each of the object frames and flat field frames. 
The bias-subtracted frames were trimmed to a 
size of 700 x 850 pixels in order to avoid the ef- 
fects due to vignetting at the edges. The pixel- 
to-pixel response variation in each of the ob- 
ject frames were corrected by dividing the object 
frames by master flat frames in each band sep- 
arately. The master flat frame in a filter band 
was obtained by averaging the best flat frames 
in the same filter and normalizing by the mean 
intensity level of the averaged frame. Cosmic ray 
events, seen as few isolated bright pixels, were 
removed by replacing them by the average inten- 
sity of four nearest neighbors. The two frames 
of NGC 4753, obtained in each band, were com- 
bined after alignment which increased the signal 
to noise ratio. The point spread functions of in- 
dividual and combined frames in each band were 
found to be similar to an accuracy better than 
0.5% which imply that the alignments were cor- 
rect. The pixel coordinates in each of the object 

5 IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of 
Universities, Inc. (AURA) under cooperative agreement 
with the National Science Foundation. The IRAF version 
2.11.1 was used. 
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image were converted into the standard equato- 
rial coordinates for the epoch 2000. The plate 
solutions were computed by fitting a quadratic 
polynomial between the known celestial coordi- 
nates and pixels coordinates of stars nearby to 
the object after projecting the celestial coordi- 
nates onto the plane tangent to the object cen- 
ter. The subsequent analysis utilized the final 
corrected and combined object images. The bias- 
corrected and flat-fielded B and R band images 
of NGC 4753 are shown in Figures | and | re- 
spectively, where the effect of extinction due to 
dust can easily be seen as the severe departure of 
the isophotes from the nearly elliptical shapes at 
the positions of dust lanes or patches. It is also 
clear that the extinction due to dust is larger in 
the B band as compared to that in the R band. 

3.1. Photometric Calibration 

Many authors (e.g. Chevalier &: Ilovaisky 1991; 
Mayya 1991; Anupama et al. 1994; Bhat et 
al. 1992) have emphasized the advantages of us- 
ing the standard star field in the "dipper aster- 
ism" region of the open cluster M 67 for photo- 
metric calibration. We have determined instru- 
mental magnitude of stars in the "dipper aster- 
ism" region of M 67 by profile fitting photometry 
using "DAOPHOT" routine within IRAF. The 
method is discussed in detail by Stetson (1987). 
To correct the instrumental magnitudes for at- 
mospheric extinction, we have made use of ex- 
tinction coefficients given by Mayya (1991), who, 
based on eight nights observation at VBT dur- 
ing January, 1991 to April, 1991, has determined 
the average extinction coefficients for the obser- 
vatory. The extinction corrected instrumental 
magnitudes and colors were transformed into the 
standard BVR system by fitting the following 
equations to the data: 

V-v = a v + (3 v (B-V) (1) 

B-V = a b - v + b - v (b - v) Q (2) 
V — R = a v - r + (3 v - r (v - r)o (3) 



B - R = a b _ r + (3 b _ r (b - r) (4) 

where the upper case letter B, V and R are the 
standard magnitudes in the corresponding filters, 
taken from Chevalier & Ilovaisky (1991). The 
linear regression model chosen to fit the data 
was ordinary least squares regression of y on x 
(OLS(Y|X)), which has been recommended for 
the data used for calibration purposes (see Isobe 
et al. 1990; Feigelson & Babu 1992). The derived 
transformation coefficients are given in Table 2. 

3.2. Surface Photometry 

We have carried out isophotal analysis of the 
galaxy NGC 4753 in B, V and R bands. The 
shapes of the isophotes were analyzed using the 
ellipse-fitting routine within the "STSDAS" f] 
software package (for details, see Jedrzejewski 1987). 
A proper background subtraction is crucial in 
the analysis. We determined the median inten- 
sity level (and its associated dispersion) in several 
25 x 25 pixels boxes in the source free regions of 
an object image. The dispersion in intensity lev- 
els in each box was found to be less than 2.4%, 
2.0%, and 2.1% in B, V, and R band images 
respectively, however, the dispersion in the me- 
dian intensity levels was found to be 0.26%, 0.4%, 
and 0.19% in B, V, and R band images respec- 
tively. We did not find any significant gradient 
in the median intensity levels in the source free 
regions across any of the object images. The sky 
background level was determined as the mean of 
the median background intensity levels in each 
band separately. The background level, thus de- 
termined, was subtracted from the object image 
in each band. The regions covered by stars and 
dust were masked and excluded from the isopho- 
tal analysis. The deviant pixels with intensity 
values 3a below or above the mean intensity level 
were clipped and the maximum accepted fraction 
of flagged pixels in an ellipse fitting was kept at 

The Space Telescope Science Data Analysis System STS- 
DAS is distributed by Space Telescope Science Institute. 
The STSDAS version 2.0.1 was used. 
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50%. The complex dust lanes and/or patches 
pass through the center of the galaxy, so it was 
not possible to determine the center very accu- 
rately. We choose the R band image, which is 
the least affected by dust-extinction, for detailed 
isophotal analysis. Starting with trial values of 
ellipticity (e), position angle (6) and ellipse cen- 
ter, an ellipse of mean intensity at a given length 
(a) of semi-major axis was fitted. The first two 
harmonics of the deviations from the trial ellipse 
were found. The best fitted ellipses were deter- 
mined after performing a minimum of 10 itera- 
tions and a maximum of 1000 iterations to min- 
imize the deviations. The third and fourth har- 
monics of the residual intensity from the best fit- 
ting ellipse were then evaluated. The procedure 
was repeated after changing the semi-major axis 
length by 10%, taking annuli and using the me- 
dian value of the pixels for sampling along the 
elliptical path. The center was then determined 
by averaging the centers of the best fitted ellipses. 
The center, thus determined, was kept fixed and 
above described ellipse fitting procedure was re- 
peated. The same center as determined for the R 
band image was also used in the isophotal anal- 
ysis of well aligned images in other bands. 

The generated radial distribution of surface 
brightness and isophotal shape parameter pro- 
files are shown in Fig. ||. The position angle 
profile reveals twisted isophotes. The ellipticity 
profile shows that the isophotes become progres- 
sively more elliptical in the outer regions. The 
parameters a3, a4, and 63, 64 are the amplitudes 
of sin39, sin^O, and cos36, cosAO coefficients, re- 
spectively, of the isophotal deviation from the 
perfect ellipticity. It should be pointed out that 
in the inner regions (a < 10" ) of NGC 4753, 
covered by complex dust lanes and/or patches, 
more than 50% of the pixels are affected by dust- 
extinction. The isophotal parameter profiles in 
the inner region, therefore, should not be taken 
very seriously. It has been tried to estimate the 
correct dust-free intensity by decreasing the ac- 
cepted fraction of pixels to fit ellipses in the inner 



regions. This may have affected the actual shape 
parameters in the central regions. The core re- 
gion of the galaxy of size of order seeing disk was 
excluded from the isophotal analysis as these are 
affected by the seeing. 

3.3. Color Index maps 

We have generated color index maps (B — R, 
B — V, V — R) using the broad band images to 
find out the distribution of dust. Instrumen- 
tal colors were corrected for atmospheric extinc- 
tion using the average extinction coefficients of 
Mayya (1991) and then converted into standard 
colors using the transformation coefficients given 
in Table 2. The B—R color index map is shown in 
Fig. superimposed on which are the contours 
of Ha image (Singh et al. 1995). The brighter re- 
gions represent the part of the galaxy which are 
redder in color and hence represent cooler /dusty 
regions within the galaxy. The B — R color in 
the dust-occupied regions within the galaxy was 
found to vary from 1.92 to 2.68. The maximum 
B — R color is found to be 2.68 at the center. In 
the dust-free regions , the B — R color as derived 
from the isophotal analysis was found to vary 
from 1.87 ± 0.07 at a semi major axis length of 
6.0" to 1.38 ± 0.07 at the semi major axis length 
of 246.7". In the dust-free regions, the color gra- 
dient with respect to logarithmic galactocentric 
radius r, , wa s found to be —0.21 ± 0.05 

for NGC 4753. This value can be compared 
with the color gradients ( ^fn^r ) —0.26, —0.23, 
-0.14 for the SO galaxies NGC 3414, NGC 3607, 
NGC 5866 respectively as derived by Vader et 
al. 1988. Hence, the average color gradient of 
NGC 4753 in the dust-free regions is normal for 
SO galaxies. Figure ||] also suggests the coexis- 
tence of dusty regions with the Ha-emitting gas. 

3.4. Extinction Maps 

Several authors have studied dust properties 
of elliptical galaxies by an indirect method in 
which dust-extinction is determined as a func- 
tion of wavelength by comparing the actual dis- 
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tribution of intensity of the galaxy with that ex- 
pected in the absence of dust-extinction. The 
dust-free intensity distribution is modeled using 
elliptical isophotes. Sahu et al. (1998) have used 
the same method to study dust properties of an 
SO galaxy NGC 2076. The projections of three 
dimensional intensity profiles of bulge and disk 
of an SO galaxy onto the plane of sky are ellip- 
tical bulge and elliptical disk isophotes. If the 
disk is inclined with respect to the plane of the 
sky, disk isophotes are more elliptical than bulge 
isophotes. In many cases, the bulge and disk 
isophotal parameters are different. In the pres- 
ence of a disk, the isophotes are not perfect el- 
lipses. The deviations of isophotes of NGC 4753 
from perfect ellipses are revealed by the presence 
of third and fourth order harmonics (o3, 63, a3, 
and 64). Therefore, the projected intensity dis- 
tribution of bulge and disk can be modeled by 
using the isophotal parameters and the harmon- 
ics. The dust-free model images of NGC 4753 
were generated by interpolating the fitted isopho- 
tal parameters (dust-free) with polynomial of or- 
der three including the third and the fourth or- 
der harmonics determined above, thus incorpo- 
rating the disk-like structure. In order to check, 
whether the dust-free intensity distribution of 
disk and bulge has been correctly modeled, we 
have carried out isophotal analysis of the dust- 
free model images of NGC 4753, and find that 
the profile of the parameter 64 derived from the 
isophotal analysis of dust-free "model image" of 
NGC 4753 is very similar to the profile of the pa- 
rameter 64 derived from the isophotal analysis of 
the "observed image" of NGC 4753. Therefore, 
we conclude that the modeled dust-free intensity 
distribution of NGC 4753 has not been affected 
significantly by the presence of a stellar disk and 
the model image represents the actual dust-free 
intensity distribution of NGC 4753. The model 
images were used to create extinction maps in 
magnitude scale as follows: 

A X = -2.5log 10 (I X ,obs/ h,model) (5) 



where A= B, V, R. Note that absolute flux cal- 
ibration is not needed for this purpose. The ex- 
tinction map of NGC 4753 in B band is shown 
in Fig. ||. The brighter areas within the galaxy 
represent regions of large optical depth associ- 
ated with the dust-extinction. The extinction in 
B band was found to vary from 0.881 ± 0.036 
to 0.438 ± 0.058. This estimate excludes cen- 
tral regions (a < 10") because of the poor accu- 
racy of the ellipse-fitting procedure. Figure |5| also 
shows overlapped contours of Ha image, which 
demonstrates that Ha-emitting regions and dust- 
occupied regions are co-spatial. The extinction 
maps in different bands were used to study the 
dust properties, as described in the next section. 

3.5. Extinction Curve 

To quantify the extinction due to dust, we se- 
lected regions in the extinction map obviously 
occupied by dust. The numerical values of ex- 
tinction, A\ (A = B, V, R), were calculated 
as the mean extinction within square boxes of 
size 5x5 pixels, which is comparable to the 
size of the seeing disk. The uncertainties associ- 
ated with the derived extinction values were esti- 
mated from the pixel-to-pixel scatter within each 
box. Since the isophotal model fits were mostly 
unreliable in the central regions (a < 10"), we 
considered safer to exclude these regions from 
the estimation of extinction values. The average 
extinction values over the dust-occupied regions 
(a > 10") within the galaxy image were assigned 
to the dust-occupied central regions. Since con- 
tribution from Ha and [Nil] line emission to the 
R band are mostly concentrated in the central 
regions (a < 10") of NGC 4753, the derived ex- 
tinction values in R band are not significantly 
affected by the presence of line emission. 

To determine the ratio of total extinction 
to selective extinction, we fitted BCES (Bivari- 
ate Correlated measurement Errors and intrinsic 
Scatter) least squares regression lines between 
different extinction values. The BCES method 
is a direct generalization of the ordinary least 
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squares (OLS) regression method, modified to 
accommodate the measurement errors either cor- 
related or uncorrelated and intrinsic scatter (for 
details, see Akritas et al. 1996). In the presence 
of errors, the use of OLS regression method can 
cause considerable bias, since it does not take 
into account measurement errors in the variables. 
Apart from the measurement errors in the de- 
rived extinction values, there are also intrinsic 
scatter depending on spatial distribution of dust 
and their properties. The measurement errors for 
extinction in two bands are uncorrelated. As a 
consequence, we must use the BCES technique 
ignoring the correlated errors in order to derive 
the linear regression lines relating extinction in 
two bands. We fitted the BCES regression lines 
A y on A x and A x on A y (x, y = B,V, R; x ^ y). 
The bisector of these two regression lines was 
chosen as the best fit line. The bisector treats 
the variables symmetrically and has been recom- 
mended for scientific problems where the goal is 
to estimate the underlying functional relation- 
ship between the variables (see Isobe et al. 1990; 
Feigelson & Babu 1992). In order to compare 
our results and to see the effects of errors on the 
regression coefficients, we have also fitted OLS 
regression lines using the method described by 
Isobe et al. (1990) between extinction values in 
different bands. The BCES and OLS regression 
coefficients are given in Table 3 and the fitted 
straight lines are shown in Figs. |6|(a) to (d). In 
Fig. ^(a), the lines marked as 1, 2, and 3 are the 
linear regression fits of Ab on Ay, Ay on ig, 
and the bisector of the two lines respectively de- 
rived using the BCES technique. In Fig. §(b), 
the lines 1, 2, and 3 are same as those in Fig. ||(a) 
but derived using the OLS method. As can be 
seen in Table 3, the lines marked 1 in Figs. ||(a) 
and (b) are significantly different. It is also seen 
in Table 3 that the lines marked 2 in Figs. ||(a) 
and (b) are significantly different but the BCES 
bisector and OLS bisector lines marked 3 in Figs. 
||(a) and (b) are similar. The difference in lines 
1 in Figs, ^(a) and (b), and also the difference 



in lines 2 can be explained as the effect of large 
errors as shown in Fig. ^(a). The BCES bisec- 
tor and OLS bisector yielded intercepts close to 
zero which are within one standard deviations as 
expected. This further suggested the validity of 
use of bisector line between extinction values in 
different bands. Henceforth, we use the BCES 
bisector line as the best fit regression line. The 
linear regression fit between An and Ab using 
the BCES bisector method is shown in Fig. ||(c). 
Similarly, the BCES bisector line of Ay and An 
is shown in Fig. |6|(d). 

The best fitting BCES bisector slopes and 
their associated uncertainties were subsequently 
used to derive the ratio of total extinction to se- 
lective extinction, R\ = Ab-A v > anc ^ their as- 
sociated uncertainties. The average extinction 
curve for the areas occupied by dust in NGC 4753 
is shown in Fig. |j] along with the extinction 
curve of our Galaxy for comparison. The to- 
tal extinction to selective extinction values for 
the Milky Way were derived in the same manner 
as for NGC 4753 from the ratios of extinction 
which, for the Galaxy, were taken from Rieke & 
Lebofsky (1985). Figure [7| shows that the ratio 
R\ = Ab-a v var i es linearly with inverse wave- 
length which is consistent with the result that 
for small grain size x < 1, Q ex t °< A -1 , where 
x = =?2, a is the grain radius and Q ex t is the 
extinction efficiency. From Fig. [7] it is obvious 
that the extinction curve for NGC 4753 is very 
similar to the Galactic extinction curve in the 
visible region. Also, the ratio of total extinc- 
tion to selective extinction, Ry = 3.1 ± 0.3 for 
NGC 4753, is similar to the Galactic value of 3.1. 
The extinction curves in Fig. [7| imply that the 
grain size distribution within the dusty regions of 
NGC 4753 is similar to that of our Galaxy. This 
conclusion about the grain size distribution can 
be used to determine the total dust mass within 
NGC 4753, as discussed below. 
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3.6. Dust Mass from Optical Extinction 

We have used the method of Goudfrooij et 
al. (1994c) to derive the dust mass from optical 
extinction values. The total dust mass can be es- 
timated by integrating the dust column density, 
given by 



ran 

J a m. t 



-7ra 3 pdn{a)da x 



(6) 



over the dust-occupied areas. In equation 6, pd is 
the grain mass density, n(a) is the grain size dis- 
tribution function, a m i n and a max are the lower 
and upper cutoffs of the grain size distribution. 
The dust column length, lj, along the line of sight 
is determined from 



A x = 1.086C«rt(A) x l d 



(7) 



where A\ is extinction in magnitude at wave- 
length A, and C ex t{\) is the extinction cross sec- 
tion of spherical grains per unit volume. C ex t(X) 
can be written as 



Ceit (A) 



Q ex t(a, \)ira 2 n(a)da (8) 



where Q e xt(a,X) is the extinction efficiency of 
grains of radii a at wavelength A. In computing 
C ex t(X), contribution of different chemical com- 
positions of the dust grains is taken into account. 

Since the extinction curve of NGC 4753 is 
quite similar to that of our Galaxy, we have em- 
ployed the grain size distribution of Mathis et 
al. (1977), i.e., 



n(a)da = Ain^a ' da 



with a r 



0.005/um and a r 



(.Q>min ^ Q> 5: O-max) 

(9) 

wx = 0.22/um in case 
Ry = 3.1. In above equation, nu is the number 
density of H nuclei, Ai is the abundance of grains 
of type i per H atom. As to the composition of 
grains, we assume spherical grains that are com- 
posed of either graphite or silicates, with equal 
abundances. This is justified by the fact that the 
Galactic extinction curve can be fitted assuming 



the above composition (Mathis et al. 1977). We 
have used parameterized V band extinction ef- 
ficiencies of Goudfrooij et al. 1994c for the two 
types of grains. The adopted values are 



Q 



ext,silicate 



and 



Q 



ext, graphite 



0.8 a/dsMcate £oi a < dsMcate 

0.8 for a > a siUcate 



2.0 o/Ograp/tite ^ ^ ^graphite 

0.8 for a > a gra phite. 



with a si i icate = 0.1/xm, and a si i icate = 0.05/um. 
The values of various parameters used for esti- 
mation of total dust mass are given in Table 4. 
The average dust mass per pixel was computed 
using equations 6 to 9, and the average value of 
Ay- We have adopted a distance of 8.7 Mpc for 
NGC 4753 (see Table 1). To find the total dust 
mass of NGC 4753, we integrated the average 
dust mass per pixel over all the dust-occupied 
pixels which were determined based on the crite- 
rion that a pixel within the galaxy image is con- 
sidered to be dusty if the pixel value in the ratio 
of observed and smooth model image is 2a or 
more below unity. The total mass of dust within 
NGC 4753, calculated thus, is 1.5 x 10 5 Mq. The 
mass of dust can also be determined based on far 
infrared (FIR) emission, as discussed in the next 
subsection. 

3.7. IRAS properties 

NGC 4753 was detected as a point source in 
the IRAS survey at all four bands 12/xm, 25pm, 
60pm and lOO^m. We obtained IRAS flux densi- 
ties from Knapp et al. (1989). The flux densities 
at 60 pm and 100 pm were corrected for the con- 
tribution of circumstellar dust emission (Goud- 
frooij et al. 1995). The dust temperature was 
determined to be 30.4 K based on the ratio of 
flux densities at 100 pm and 60 pm under the as- 
sumption that the FIR emission from the dust 
grains within NGC 4753 is governed by an emis- 
sivity law where emissivity is proportional to A -1 
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at wavelengths 200/um (Schwartz 1982; Hilde- 
brand 1983; Kwan & Xie 1992). Since a distri- 
bution of temperature may be more appropriate 
for dust, the derived dust temperature should be 
regarded as a representative value. We have esti- 
mated the dust mass of NGC 4753 following the 
method outlined by Hildebrand (1983), Thron- 
son et al. (1986) and Goudfrooij et al. (1995) 
which is based on measurements of FIR flux den- 
sity /(;/), dust temperature T^, dust emissivity 
Q{f) and size of the grains. The total dust mass 
is given by 



M d 



B{v,T d )'Q(y 



Pd 



(10) 



where p d is the specific grain mass density, a is 
the average grain radius weighted by the grain 
volume and D is the distance of the galaxy. The 
derived cool dust mass is 3.46 x 1O 5 M . As- 
suming that the infrared spectrum of NGC 4753 
can be fitted by a function of the form f{y) oc 
vB(v,T<i) which is a good approximation to the 
total infrared spectrum of most galaxies (Telesco 
& Harper 1980; Rickard & Harvey 1984), the to- 
tal infrared luminosity can be written as (Thron- 
son et al. 1986) 

Lir(L q ) = 3.7xW- 12 D 2 f(v)T%\ 4 [exp(lA4xW 4 /\T)-l} 

where D is the object distance in Mpc, the flux 
density f(u) is in Jy, and A is in /im. Using 
the above relation the total infrared luminosity 



of NGC 4753 is estimated to be 6.6 x 10 8 L 



©■ 



4. Discussion 



The surface brightness profile and isophotal 
shape parameters (Fig. ||) of NGC 4753 reveal 
some intensity features of this early-type galaxy. 
The parameter 64, the amplitude of cos4# coef- 
ficient of the isophotal deviation from the best 
fitting ellipse, is interpreted as being due to an 
embedded disk (64 > 0) or due to boxiness of 
isophotes (64 < 0) (e.g. Peletier et al. 1990). 
The positive values of 64 indicate a disk compo- 
nent in NGC 4753. The galaxy is very flattened 



in its outer parts. The presence of disk compo- 
nent is supported by large rotational velocities, 
v > 250km s _1 measured by Chromey (1973). 
The disk component and the twisted isophotes as 
revealed by the position angle profile, Fig. |3|(b), 
are supported by a model where an inclined disk 
twisted by differential precession fits the observed 
dust lanes lying on the disk very well (Steiman- 
Cameron et al. 1992). 

We have derived the dust mass of NGC 4753 
from optical extinction as well as from IRAS flux 
densities. The ratio of the two dust masses, 
jgjf^f; is 2.28. The discrepancy can be at- 
tributed to the fact that the dust mass derived 
from optical extinction is a lower limit because 
the extinction is observed only when the dust is 
in front of the stars. The dust is probably dis- 
tributed between the stars or stellar system, the 
actual dust mass is expected to be higher than 
the dust mass derived from optical extinction. 
Goudfrooij et al. (1995) have reported that the 
average ratio / ^ d > IRAS ) to be 8.4±1.3 for a large 

° N Md, optical ' ° 

sample of elliptical galaxies for which the pres- 
ence of dust is revealed by both FIR emission 
and optical dust lanes or patches. This ratio is 
significantly higher than the ratio for NGC 4753 
which implies that the dust mass discrepancy in 
NGC 4753 is not as severe as for an average el- 
liptical galaxy. Sahu et al. (1998) have reported 



that 



Ma. 



1.14 for an SO galaxy NGC 2076. 



M d ^optical 

This indicates that the difference in dust mass 
derived from FIR emission and optical extinc- 
tion is not significant for SO galaxies in contrast 
to the dust mass discrepancy in elliptical galax- 
ies. This is an expected result because the dust 
masses derived for spiral galaxies from optical 
and IRAS data reveal that M^iras ^ ^d,optical 
(Goudfrooij 1996) and the SO galaxies, being in- 
termediate in the morphological sequence from 
ellipticals to spirals, can be expected to have in- 
termediate properties. 

A rough estimate of molecular gas content of 
NGC 4753 can be made from the dust mass de- 
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rived from FIR flux if we assume that the ratio 
of molecular gas to dust mass to be in the range 
50 - 700 as given by Wiklind et al. (1995) for 
emission-line regions . This would render the 
molecular gas mass of NGC 4753 to be in the 
range 0.2- 2.4 x 1O 8 M . The masses of different 
contents of NGC 4753 are given in Table 5. 

The origin of the observed interstellar dust can 
be either internal (e.g. due to mass loss from 
stars) or external (galaxy-galaxy interaction or 
merger). In the solar neighborhood, the potential 
sources of dust grains are supernovae, red giant 
stars, protostellar nebulae and planetary nebu- 
lae with relative weights of 100, 37, < 16, and 
6 respectively (Dwek Sz Scalo 1980). In the case 
of early-type galaxies, the rate of star formation 
and hence the rate of supernova explosions are 
much lower than the corresponding rates in spi- 
ral galaxies. The rate of supernova explosions in 
early-type galaxies is about a factor of 16 smaller 
then the rate in spiral galaxies (Sbc-Sd) (Cappel- 
laro et al. 1993). Therefore, it can be expected 
that the dominant sources of dust grains in early- 
type galaxies are the red giant stars. In the fol- 
lowing, we consider dust mass within NGC 4753 
arising from red giant stars and neglect the other 
possible sources of dust grains. The total mass 
loss rate from red giant stars within NGC 4753 
can be estimated from the excess infrared (IR) lu- 
minosity arising from circumstellar shells of the 
red giant stars (see Soifer et al. 1986). The excess 
IR luminosity arising from circumstellar shells of 
red giant stars can be determined by comparing 
the IR luminosities from photospheres of normal 
stars and red giant stars and by counting the to- 
tal number of red giant stars within NGC 4753. 
Since individual stars within NGC 4753 could not 
be resolved, it is not possible to determine the 
excess IR luminosity arising from circumstellar 
shells of the red giant stars. Hence the total mass 
loss rate from red giant stars within NGC 4753 
cannot be determined directly. However, follow- 
ing the analysis of Jura et al. (1987) and assum- 
ing that stellar population within NGC 4753 is 



similar to that of elliptical galaxies or bulge of 
M 31, the total mass loss rate from red giant 
stars within NGC 4753 can be estimated. The 
total mass loss rate from red giant stars within 
an elliptical galaxy can be written as (Jura et 
al. 1987) 

M (M yr- 1 ) = 4.3xl0~ 30 ^(12/xm) (erg s^Hz' 1 ) 

(11) 

where Lj,(12/im) is the luminosity per unit fre- 
quency at 12/im. For NGC 4753, we calculate 
L„(12/im) to be 3.0 x 10 26 erg s" 1 Hz' 1 from 
the measured flux at 12/im and adopting a dis- 
tance of 8.7 Mpc (see Table 1). Substituting 
L 1/ (12/im) into above equation, we estimate the 
total mass loss rate from red giant stars within 
NGC 4753 to be 1.3xlO- 3 M yr' 1 . The derived 
mass loss rate is uncertain because stellar pop- 
ulation within elliptical galaxies and NGC 4753 
may be significantly different, however, the mass 
loss rate is not expected to vary by large fac- 
tors. If we assume that the gas-to-dust ratio in 
the circumstellar shells of red giant stars to be 
100 : 1 by mass, with a factor of 2 — 3 uncer- 
tainty (see Knapp et al. 1993), the total loss rate 
of dust mass by red giant stars within NGC 4753 
is about 1.3 x 10~ 5 Mq yr^ 1 . The mass of dust 
accumulated within NGC 4753 can be estimated 
by the mass loss from stars and by the mech- 
anisms of destruction of dust (see below). The 
rate at which dust mass is accumulated within 
NGC 4753 can be written as 

Z-MMA-mm^ (i2) 

where d ^'" is the loss rate of mass in the form 
of dust by red giant stars and r _1 is the de- 
struction rate of dust. The second term in the 
right represents the rate at which mass of dust 
within NGC 4753 is destroyed. Assuming that 
the mass loss from red giant stars begins at an 
age ~ 10 6 yr of NGC 4753 and taking r _1 ~ 
7.14 x 10~ 10 yr^ 1 (see below), we have solved nu- 
merically the above equation for the mass of dust 
accumulated as a function of time. The build up 
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of mass of dust within NGC 4753 is shown in Fig. 
||. Assuming an age of 10 10 yr for NGC 4753, the 
total mass of dust accumulated within NGC 4753 
is calculated to be 1.8 x 10 4 Mq, which is about 
a factor of 10 lower than the measured mass of 
dust. This indicates that the mass loss rate from 
red giant stars is probably not sufficient to ac- 
count for the observed dust within NGC 4753. 
Based on our simplified calculation given above, 
the possibility of an external origin of dust within 
NGC 4753 can not be ruled out. 

X-ray emission of varying amounts is seen 
in several early-type galaxies. Integrated X-ray 
emission from discrete sources and the presence 
of hot gas (~ 10 7 K), both can contribute to the 
overall X-ray emission. Hot plasma, if present, 
can destroy the dust very efficiently (Draine & 
Salpeter 1979a). The ratio of total X-ray lumi- 
nosity to absolute blue luminosity is a good indi- 
cator of the presence of hot gas in an early-type 
galaxy (Canizares et al. 1987). Therefore, in or- 
der to estimate the dust destruction time scale, 
we examine X-ray emission from NGC 4753. The 
quantity log(j^) (Lx in ergs s _1 and Lb in Lq) 
for NGC 4753 is 29.42, based on which Canizares 
et al. (1987) would conclude that the total X- 
ray luminosity of NGC 4753 is consistent with 
that expected from the integrated X-ray lumi- 
nosity of discrete sources. The galaxy NGC 4753 
is a member of the lowest group in Kim et 
al. (1992) where it has been shown that early- 
type galaxies in the lowest group have a very 
soft X-ray excess which amounts to about half 
the total X-ray emission. This emission could 
be due either to a cooler ISM or to the inte- 
grated emission of soft stellar sources. In any 
case, ISM of these galaxies is not heated to high 
X-ray temperatures. This implies that the dust 
in NGC 4753 is unlikely to be embedded in a 
very hot plasma. In the absence of hot plasma, 
the most effective destruction mechanism for re- 
fractive grains (e.g., graphite and silicate) are 
sputtering and grain-grain collisions in low veloc- 
ity shocks {v shock < 60 km s~ l ), and sputtering 



in supernova driven blast waves (Goudfrooij et 
al. 1994c). Using the model of Draine & Salpeter 
(1979b), appropriate for X-ray faint early-type 
galaxies, and recent estimate for the rate of su- 
pernova explosions in early- type galaxies, Goud- 
frooij et al. (1994c) have estimated the life time 
1.4 x 10 9 yr for 0.1 //m refractory grains. We 
have used this value to calculate the mass of dust 
accumulated within NGC 4753, given above. 

Ha+[NII] emission line regions in the center 
of NGC 4753 are embedded in the regions with 
substantial optical extinction (see Figs. ^| and 
||).This association of emission line regions and 
regions with substantial dust absorption points 
towards the coexistence of dust and warm gas 
(~ 10 4 K) within NGC 4753. The observed 
Ha line emission can be accounted by photo- 
ionization from post asymptotic giant branch 
stars (see Singh et al. 1994; 1995), which can 
also supply significant amounts of dust. 

As to the fate of cool dust and gas, we ex- 
amine the possibility of star formation within 
NGC 4753. The galaxy occupies an intermedi- 
ate position but close to the lower right end in 
the phenomenological IRAS color-color diagram 
of Helou (1986), where both the cirrus compo- 
nent and active star forming regions contribute 
to the total FIR emission. However, the FIR 
emission from active star forming regions is less 
than 50% of the total FIR luminosity, so that the 
total FIR emission of NGC 4753 can not be in- 
terpreted as due to current star formation. Fol- 
lowing Thronson & Telesco (1980), the current 
star formation rate averaged over past 2 x 10 6 yr 
can be calculated using the relation 

M FIR = 6.5 x 10- 10 L FIR (L Q ) (13) 

where Lpm is the total FIR luminosity. How- 
ever, considering the phenomenological model of 
Helou (1986), only the FIR luminosity due to 
the active star forming regions must be used in 
above expression. We estimate the current star 
formation rate of NGC 4753, averaged over past 
2 x 10 6 yr to be less than O.llMQyr' 1 . As dis- 
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cussed above, the cirrus component of dust con- 
tributes more than 50% to the total FIR lumi- 
nosity, hence significant amount of dust within 
NGC 4753 is in the form of cirrus. This sug- 
gests that the cirrus clouds are not destroyed ef- 
ficiently because of the lack of hot ISM. Also, 
these regions are not sufficiently dense to lead to 
gravitational instability to form stars. 

5. Conclusions 

Our isophotal analysis of NGC 4753 has re- 
vealed a twisted disk component in it. The dust 
grain properties of NGC 4753 appear to be simi- 
lar to that of our Galaxy based on comparison of 
the extinction curves. The derived mass of cold 
dust in NGC 4753 is 1.5 x 1O 5 M of NGC 4753 
based on optical extinction, and 3.46 x 10 5 Mq 
based on FIR fluxes. The discrepancy between 
the two dust masses for NGC 4753 is not as se- 
vere as the discrepancy generally found for the el- 
liptical galaxies. The accumulated mass of dust 
within NGC 4753 from mass loss by red giant 
stars after taking into account the efficient de- 
struction processes is ~ 1.8 x 10 4 Mq, which is 
about a factor of 10 lower than the measured 
dust mass. Dust and ionized gas coexist within 
NGC 4753. The current star formation rate of 
NGC 4753, averaged over past 2 x 10 6 yr, is de- 
rived to be less than 0.21 Mq^t- -1 using the FIR 
emission arising from active star forming regions. 
A significant amount of dust within NGC 4753 
seems to exist in the form of cirrus. 
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Fig. 1. — Bias-subtracted, flat-fielded and cos- 
mic ray removed B band image of NGC 4753 
superimposed on which are the contours of same 
image. The contour levels are drawn at 4.2%, 
4.5%, 5%, 6%, 7%, and 8% of the peak inten- 
sity. The outermost contour is drawn at a sur- 
face brighness 23.1 magnitude per arcsec 2 . The 
brighter lanes or patches represent dust-occupied 
regions. 

Fig. 2. — R band image of NGC 4753 overlapped 
with the contours of the same image. The con- 
tour levels are drawn at 3%, 3.3%, 4%, 5%, 6%, 
7% and 8% of the peak intensity. The outer- 
most contour corresponds to the surface bright- 
ness 22.1 magnitude per arcsec 2 . 



Fig. 3. — Results of R band isophotal shape 
analysis of NGC 4753. (a) Surface brightness 
profile, (b) position angle profile, (c) ellipticity 
profile, (d) amplitude, a3, of residual sinW co- 
efficient of the isophotal deviation from perfect 
ellipse as a function of semi-major axis length 
(SMA). (e) amplitude, a4, of residual sinAB co- 
efficient of the isophotal deviation from perfect 
ellipse as a function of semi- major axis length 
(SMA). (f) amplitude, 63, of residual cos39 co- 
efficient of the isophotal deviation from perfect 
ellipse function f semi-major axis length 

(SMA). (g) amplitude, 64, of residual cosAO co- 
efficient of the isophotal deviation from perfect 
ellipse as a function of semi- major axis length 
(SMA). 



This 2-column preprint was prepared with the AAS IATj^X 
macros v4.0. 



Fig. 4. — B — R color image of NGC 4753 super- 
imposed on which are the contours of Ha image. 
The brighter regions represent part of the galaxy 
which are redder in color. The maximum B — R 
color is found to be 2.677 at the center. The 
Ha+[NII] contour are drawn at levels 5%, 7.5%, 
10%, 15%, 25%, 50% and 80% of the peak inten- 
sity. The Ha+[NII] image is taken from Singh et 
al. (1995). 
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Fig. 5.— Extinction image of NGC 4753 
in B band. The brighter areas represent re- 
gions of large optical depth associated with the 
dust-extinction. The superimposed contours of 
Ha+fNII] (Singh et al. 1995) map are drawn at 
levels 5%, 7.5%, 10%, 15%, 25%, 50% and 80% 
of the peak intensity and show spatial association 
of dust and Ha gas. 



Fig. 6. — BCES and OLS regression lines be- 
tween extinction values in different pass bands, 
(a) BCES regression lines between extinction in 
V and B bands. The line 1 is the BCES regres- 
sion line of Ay on ig, the line 2 is the BCES 
regression line of Ab on Ay, and the line 3 rep- 
resents the BCES bisector of the above two lines 

1 and 2. (b) OLS regression lines between Ay on 
Ab- The line 1 is the OLS line of Ay on Ab, line 

2 is the OLS line of Ab on Ay, and line 3 is the 
OLS bisector of lines 1 and 2. (c) BCES bisector 
regression line of Ab and Ay. (c) BCES bisector 
regression line of Ar and Ab- The quantities A\ 
(A = B, V, and R) represent the dust-extinction 
in the band A. 



Fig. 7.— Extinction curve of NGC 4753. The 
Galactic extinction curve is also plotted for com- 
parison. The quantity R\ (= aJ^a v ) is the ratio 
of total extinction to selective extinction in the 
band A = B, V, R. 



Fig. 8. — The total dust mass accumulated 
within NGC 4753 as a function of time. 



15 



